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Abstract 



^ ' We study a Penrose-Fife phase transition model coupled with homogeneous Neumann bound- 

ary conditions. Improving previous results, we show that the initial value problem for this model 
admits a unique solution under weak conditions on the initial data. Moreover, we prove asymptotic 
regularization properties of weak solutions. 
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1 Introduction 



5h 

The Penrose-Fife system, proposed by O. Penrose and P. Fife in 26, 27 , represents a thermodynam- 
ically consistent model for the description of the kinetics of phase transition and phase separation 
processes in binary materials. It couples the singular heat equation (|2.3p for the absolute temperature 
i? with a nonlinear relation describing the evolution of the phase variable X which represents the local 
proportion of one of the two components. This can be of the fourth order in space (cf. (I2.6[) - (|2.7[) 
below), in case the physical process preserves the total mass of X (conserved Penrose-Fife model, de- 
scribing phase separation) or of the second order in space (cf. (I2.5P below) , in case the total mass of X 
is admitted to vary (non-conserved Penrose-Fife model, describing phase transition). In the conserved 
case, the equation for X is usually written as a system by introducing an auxiliary variable w called 
chemical potential. We refer to the next section for a detailed presentation of the equation and to the 
papers jTQl [TTJ [33] for further mathematical background. 

Due to physical considerations, it is generally accepted to consider no-flux boundary condi- 
tions for X and, in the conserved case, also for w. On the other hand, various types of boundary 
conditions (for instance, no-flux, non-homogeneous Dirichlet, or Robin conditions) make sense for 
the heat equation (|2.3p . which give rise to different mathematical scenarios. We refer the reader to 
[TTJ U2 HI 121 US [33] for the case of Robin (or "third type") conditions, to [HQ13 for the Dirichlet 
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case, and to [TU] [HI HE] for the homogeneous Neumann case, which is probably the most difficult one 
due to lower coercivity properties of the elliptic operator in (|2.3I) . 

Our aim in this paper is that of improving existing results on the homogeneous Neumann 
problem for the Penrose-Fife model both in the non-conserved and in the conserved case. Actually, 
our results will cover both situations, generally with minor variations in the proofs. As a first property, 
we will show that the problem admits a unique solution under weak assumptions on the initial data. 
Actually, noting that the system admits a natural Liapounov functional representing the total energy, 
we will prove that a (unique) weak solution exists for any initial data ($o, ^o) having finite energy. This 
improves existing results which assume some extra summability condition on $o, typically $q £ L 2 (Q). 
As a side effect, we pay the price that the heat equation has to be interpreted in the generalized (H 1 )'- 
framework developed by Damlamian and Kenmochi in [13j . Namely, the (thermal part of the) energy 
has to be intended as a (relaxed) functional operating on the negative order Sobolev space if 1 (0)' 
(cf. (|2.20[) below) and also the relation linking -d to its inverse has to be stated properly. On the other 
hand, if we know in addition that $o £ L 1 (fi), then we can prove that i?(i) £ L x (f2) for all t > 0; 
moreover, both the energy functional and the relation between $ and its inverse can be written in the 
usual (pointwise) sense (cf. (|2.19|) below). 

In the subsequent part of the paper, we prove our main results, which regard uniform time- 
regularization properties of weak solutions. In this frame, we will actually present two theorems. 
Firstly, we will show that, for any T > 0, there exists a constant $ > depending only on the 
"energy" of the initial data and on T such that d(x,t) > $ for a.e. (x,t) £ (T, +oo) x Q, Second, 
we will prove that a similar bound from above (i.e., i!)(x,t) < ■& for a.e. (x,t) £ (T, +oo) x fi) and a 
suitable > 0) holds provided that the initial temperature $o 5 in addition to the "energy" regularity, 
satisfies the additional hypothesis $o £ L 3+£ (il) for some e > 0. This additional condition appears 
also in other works concerning L°°-regularization properties of the solutions to very fast diffusion 
equations like (12.31) on M 3 (see, e.g., [33], [5] and the references therein). In particular, in three space 
dimensions, the exponent p = 3 happens to be critical for the boundedness (for strictly positive times) 
of the solutions to (I2.3[) : starting from initial data in L P (M. 3 ), p > 3, implies boundedness of the 
solutions for strictly positive times (see [5]). For p < 3, the situation is drastically different, as the 
self similar solution (|2.40[) shows. 

We finally note that the proved uniform bounds permit, by standard methods, to improve fur- 
therly the regularity of solutions for strictly positive times. In particular, our estimates complement 
a recent paper by Priiss and Wilke [28 who show maximal regularity estimates for the conserved 
Penrose-Fife model under the conditional (i.e., unproved in their paper) assumption that the temper- 
ature ■& satisfies the uniform bounds $ < $(x, t) < -d almost everywhere. Thanks to our results, the 
maximal regularity estimates of Priiss and Wilke hold for all strictly positive times and all weak solu- 
tions emanating from initial data satisfying the "energy regularity" plus the condition do £ L 3+e (Sl). 
It is also worth noting that the uniform bound -d > $ implies, in the nonconserved case, the so-called 
"separation from singularities" property for X in case the configuration potential of the system (i.e., 
the function b defined in (I2.10[l ) has a bounded domain. For instance, in the physically relevant case 
of the logarithmic potential 

b{r) = (1 + r) log(l + r) + (1 - r) log(l - r), (1.1) 

whose domain is [—1, 1], this means that — 1 + S < X(x,t) < 1 — S for some 8 > and for all times 
t > 1. Unfortunately, we cannot prove such a property in the case of the conserved model. Actually, 
up to our knowledge, this is unknown, at least in the three-dimensional case, also in for the (simpler) 
conserved phase- field model of Caginalp type (cf. the related discussion in [13] )• 

The remainder of the paper is organized as follows. In the next section we will present our 
hypotheses and state our results. The proofs will be detailed in the subsequent Section [3] 

Acknowledgment. The authors would like to thank Prof. Philippe Laurengot for discussions re- 
garding the strategy of some proofs. 
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2 Main results 

Let fl be a smooth bounded domain of M. d , d 6 {2, 3}. For the sake of simplicity, let us assume |f2| = 1 
so that ||u[[jy>(n) < ||«||i9(n) for all 1 < p < g < +oo, v s L q (Q). For simplicity, we will often write 
|| • ||p in place of || • ||z,P(n)- Let H := L 2 (Q), endowed with the standard scalar product (•, •) and norm 
| • ||. Let also V := H 1 ^). We note by || • \\x the norm in the generic Banach space X. and by (•, -)x 
the duality between X' and X. 

For any function, or functional z, defined on O, we can then set 

=" : =H L-ls (21) 

where the integral is substituted with the duality (z, 1) in case, e.g., z <E V . 
We also define the elliptic operator 

A:V->V, (Av,z) := I VvVz. (2.2) 

Then, for t S (0, +oo), we consider the singular heat equation 

§ t + Au = -X u in V, (2.3) 
u = — — , almost everywhere in 57. (2-4) 

In the non- conserved case, this is coupled with the equation 

X t + AX + b(X) - X = u, in H, (2.5) 

while in the conserved case, (|2.3D - (|2.4[) is coupled with the system 

X t + Aw = 0, in V', (2.6) 
w = AX + b(X) - X - u, in V . (2.7) 

In both cases, we will take the initial conditions 

0|t=o = i?o, X\t=o = X , in n. (2.8) 
The nonlinear function b is assumed to satisfy 

b 6 C°'Ul;R), 6(0) =0, b'(r) > a.e. in I, lim 6(r) signr = lim b'(r) = +oo, (2.9) 

where /, the domain of b, is an open, possibly bounded, interval of K containing 0. We also set 



b(r) := / 6(a) ds (2.10) 



o 

in such a way that 6 is a convex function satisfying 6(0) = 0. 

Concerning the initial datum Xq, we will always assume 

X eV, 6(X )eL 1 (fi). (2.11) 

To specify the regularity required for the initial temperature $o, we first need to introduce some 
convex analysis machinery. First of all, we set 

j(«):=-log(-«), f(z):=-l-logz, (2.12) 

respectively for v < and z > 0. Then, j and j* can be seen as a couple of conjugate functions 
according to the standard theory (cf., e.g., [7 ). This permits to introduce the convex functional 

J:V -> [0,+oo], J(«) := / (2.13) 
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where we have implicitly set J(v) = +00 for those v G V such that j(v) is not summable (this 
happens, for instance, when v is strictly positive on a set of strictly positive measure). The set 
{v £ V : J(v) < +00} is called domain of J. Then, the conjugate of J is given by 

J* : V' (-00, +00], J*(C) := sup { (C, v) — J{v), v eV}. (2.14) 

We can also introduce the subdifferential of J with respect to the duality between V and V. Namely, 
given v G V, we set 

C€dv, v >J(v)<=^(C,z-v)<J(z)-J(v) VzeV. (2.15) 

In general, the "weak" subdifferential dyyJ is a multivalued maximal-monotone operator. Its struc- 
ture is analyzed in several papers (see, e.g., [S]; see also [5] and [TB] for the slightly different situation 
where V is substituted by Hq(Q)). Actually, it in not difficult to prove that 

(ei 1 ^), C = f(v)= a.e. in ft => ( G dv,v>J{v). (2.16) 

On the other hand, for given v G V, the set dv,yJ(v) needs not be contained in L 1 (ft). More precisely, 
a generic element £ G dyyJ{y) is a bounded measure which may have a singular part Q s . A precise 
characterization is given in [8j Thm. 3] for the -ffg-case. However, it is easy to realize that such a 
characterization extends to the case of V = i? 1 (ft) at least provided ft is smooth, the only difference 
being that Cs can be supported also on the boundary 9ft. On the other hand, if we are able to prove 
that the singular part of some element C is 0, then we still have pointwise inclusion. Namely, there 
holds that 

CGd V V'J(v)nL 1 (Sl)=>( = — a.e. in ft. (2.17) 

v 

That said, the minimal regularity required on i?q is given by 

4 e V", J*(tf )<+oo. (2.18) 
Then, we can define the "strong" energy functional of the system as 

£ s {d,X) :=^(tf + X-l-logtf+i|VX| 2 +6(X)-ix 2 ), (2.19) 

However, since in general $0 is just an element of V , it does not make sense to compute £ s ($o,Xo). 
Actually, we have to relax £ s , by defining 

£{d, X) := (i?, 1) + J*(0) +f n (*+ ||VX| 2 + b(X) - \* 2 )- (2.20) 

Actually, assumptions (I2.18[) and (12. lip are equivalent to asking that the "relaxed" energy Eo := 8(0) 
is finite. In the sequel we shall often write £(t) in place of £(d(t),X(t)). Notice that the last (I2.9[) 
ensures the coercivity of £. Notice also that the no-flux conditions entail some conservation properties. 
Actually, for the conserved system, testing (|2.3[) and (|2.6p by 1, we immediately get 

0n(*) = (0o)n, Xn(t) = (Xo)n, Vt > 0, (2.21) 
while of course for the non-conserved system we only have 

0n(t)+Xn(t) = O?o)n + (X o )n, Vi - ' (2 ' 22) 

The following result, stating existence of strong solutions, is well-known (see, for instance, [T5] for a 
proof): 

Theorem 2.1. Let us assume (123)) . f2~Tg]t . (|2TTT|) and, additionally, let 

tf € H. (2.23) 
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Then, in the non-conserved case, there exists a unique triplet u, X) satisfying, [or all T > 0, 

tf-logtfe-L 00 ^,! 1 ;!, 1 ^)), e L 2 {0,T;H), (2.24) 
ueL 2 (0,T;F), (2.25) 
X e H 1 (0,T;H)nL oo (0,T;V)DL 2 (Q,T;H 2 (n)), b(X) E L 2 (0, T; H), (2.26) 

and solving J2.3|) - (|2.4|) and (|2.5p a.e. in (0, T), together with the initial conditions (|2.8|) . Ana7ogous7y, 
in the non-conserved case, we have a unique triplet ($,w, X) satisfying, for all T > 0, (|2.24p - (|2.25[) . 
together with 

XeH 1 {Q,T;V')r\L co {Q,T;V)r\L 2 {Q,T;H 2 {VL)), b(X) e L 2 {0,T: H), (2.27) 
solving dH31)-(IlIU) and dMD-(E3D a.e. in (0,T), and satisfying (gU) . 
Remark 2.2. It is worth noting (cf. [5] for more details) that (|2.23[) entails 



J*(z?o) = / (-1-log^o). (2.28) 
Thus, in fact no functional on V' appears in the above "strong" formulation. 

In the next result we will provide an existence theorem working without the additional regularity 
(|2.23p . We will pay the price of the occurrence of relaxed functionals. Moreover, we will also need to 
intend relation (|2.4[) in the relaxed sense of |13| . 

Theorem 2.3 ((V)- solutions )- Let (EH), (|2~T8| and ([2~TT|) hold. Then, there exists a unique triplet 
(•d,u,X), satisfying, for all T > 0, 

e H 1 (Q, T; V'), J*(§) S L°°(0,T), (2.29) 

together with (|2.25p and, a.e. in (0, T), equation (|2.3p . Moreover, in the non-conserved case, (|2.26p and, 
a.e. in (0, T), equation (|2.5I) hold, while, in the conserved case, (|2.27p and, a.e. in (0, T), equations (|2.6I) - 
(|2.7p hold. Moreover, there hold the initial condition (12. 8p and, a.e. in (0,T), the weak identihcation 
property 

tied v ,v>J(u). (2.30) 
Finally, any V' -solutions satisfies, for all t g [0, T], the energy equality 

4 (iiv^ii 2 + nx t |i 2 ) = j*(tf ) - j*(i?(t)) + (i? , i) - mi i) 

(Xo + ^IVXoI 2 + KXo) - i^) - J (X(t) + \\VX{t)\ 2 +b(X(t)) - \x{t) 2 ) , (2.31) 

in the non-conserved case. In the conserved case, the same holds provided that the term \\Xt\\ 2 on the 
left hand side is replaced by ||Vw|| 2 . 

We now turn to discussing regularization properties of solutions. The first result regards the function 
•u. 

Theorem 2.4. Let (|2"T9"]). (|2~T51) and (12~TT]) hold. Then, the V -solution either to the non-conserved 
or to the conserved problem satisfies 

IMU~(l,oo;V) + IMlL~((2,oc)xQ) < Q(E ), (2.32) 

||X]U»(i,oo ; ^(n)) + ||KX)IU~(i,oo;H) < Q(Eo). (2.33) 

Moreover, in the non-conserved case we aiso have the "separation property" 

\\b(X)\\ L oo as>oo)xn) <Q(E ). (2.34) 

Here and below, Q is a computable nonnegative function, monotone in each of its arguments, whose 
expression is independent of initial data and of time. 
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Remark 2.5. It is worth noting that the above regularization properties are in fact instantaneous. 
Indeed, with minor modification in the proofs one could easily see that (|2.32l) - (|2.33l) and (|2.34l) hold 
starting from any r > (and not only from r = 1 or 2 or 3). Of course, then the functions Q on the 
right hand sides would also monotonically depend on r _1 and possibly explode for r \0. The same 
considerations hold also for what is proved in Theorem 12.71 below . 

In case the initial temperature, beyond satisfying (|2.18[) . is an ./^-function, we can say something 
more precise on regularity: 

Theorem 2.6 ((V n i 1 )-solutions). Let f2l)j) . (j2~18l) and (|2~TT1) hold. Let also 

tfoEL 1 ^). (2.35) 
Then, the triplet u, X) given by Theorem 12.31 additionally satisfies, for all T > 0, 

e C°{[0,T];L l {tt)) (2.36) 
and the strong identification property (12.41) . 



It is easy to show that, if do € L p (fl) for p > 1, then §(t) remains in L P {VL) for t > 0. What is more 
interesting is that, if p > 3, then ■& is asymptotically uniformly bounded: 

Theorem 2.7. Let ([2T5]) . ([2TTg)l ajid (|2"TIj) hoW and Jet afeo 

i?o € L 3+£ (n) for some e > 0. (2.37) 

Then, the (V (~l L 1 )-solution either to the non-conserved or to the conserved problem satisfies the 
additional bound 

Wi- ((2 , oo)xn) <g(E 0j ||i?o||a+ e ). (2-38) 

Remark 2.8. Relation (|2.37l) suggests that 3 should play the role of a critical exponent for equation 
(|2.3p in space dimension 3. Actually, it is easy to check that the related "very-fast diffusion" equation 



i9 t + At? _1 =0 (2.39) 
over the whole space (0, +oo) x M 3 admits the similarity solution (see, e.g., [34] ) 

tf(t,x) = 2(T " f) + 2 , (2.40) 
\x\ 

which belongs to Lf oc (R 3 ) for allp < 3 and alH > and does not exhibit any instantaneous regularizing 
effect (of course, it satisfies a delayed regularization property since it extinguishes in a finite time; 
however, this effect is not expected to hold in the case of a finite domain when we have conservation 
of mass). However, we do not know what happens in the critical case of an initial datum -d belonging 
to L 3 {n). 

3 Proofs 

All proofs will be in principle given only for the conserved case which is, actually, more difficult. The 
properties holding only for the non-conserved case (as well any significant differences in the proofs) 
will be remarked on occurrence. 

3.1 Proof of Theorem 1231 

We start by proving existence for weak initial data. Given f?o satisfying (|2.18p . we then set, for n G N, 

O) n+n- 1 A0 o ,n = 0o. (3.1) 
The properties of this approximation deserve to be stated in a lemma. 
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Lemma 3.1. Let (|2.18p hold and let i? ,n be defined by (f3~Tj) . Then, i\n G V for all n G N. Moreover, 
#o,n -> $o strongfy in 7' and J*(i? 0) n) < J*{$o) Vn G N. (3.2) 
Moreover, if a&o (|2.35|) holds, then we also have 

do,n^$o strongly in L 1 (SI) (3.3) 

and 

log^o,„ eL\n), /-l-logi? ,n< /-l-logtfo, VneN. (3.4) 

Proof. Being $o G V, it is clear that i^o.n £ ^ for all n. Moreover, the ^'-strong convergence 
in (|3.2p can be proved by standard Hilbert techniques. In order to complete the proof of (|3.2[) . we 
introduce, for any given n G N, the strictly positive sequence $o,n,fc '■= max{$ ,nj j}, for fc G N. 
Correspondingly, we set wo, n ,fc := — ^ 1 - ■ Note that, by construction, Uo,n,fc G V and the map 
$o.n l— ► $o,n.k is monotone. As a consequence, we can write 

> ^ - -Al?o,»», "0,n,fc^> = (l?0,n - $0, "0,n,fc) = ($0,n,fe ~ $0, Mo.n,*) + (l?0,n ~ #0,n,fe) "0,n,k}- (3-5) 

Now, since $o,n,fc = ~ M[) 1 fc almost everywhere in f2, we have $o,n,fc G cH^y J(uo,n,fc)- Equivalently, 
"o.n.fe G <9y,y J*($o,n,fc) (where the subdifferential acts now in the duality between V and V) for any 
fc. Hence, 

(l?0,n,fc - ^0,«0,n,fc> > J*0V,fc) - J*(l?o), (3.6) 

by definition of subdifferential. On the other hand, since i9o,n,fc = i?o,n in f2 R {$o,n > l/^}> 



($0,n - #0,rj,fc,"0,ra,fc) = / ' (#0,n ~ $0,n,fc) > 0. (3.7) 

Jnn{i9 ,r.<i/fc} 
Thus, collecting the above computations we have 

J*($0,n,fe) < J*(i?o). (3-8) 

Finally, since for any n G N we have that #o,n,fc fc ^ +00 > $ 0n strongly in L p (fl) for any p G [1,6) 
(hence, a fortiori, in V), we have the following chain of inequalities: 

J*0?o,») < liminf J*(0 Ol n,*) < limsu P ■H'W) < J*(#o), for any n G N, (3.9) 

fe/M-oo fe/'+oo 

i.e., d321) holds. 

Now, we assume that also (|2.35l) holds, namely we assume that tDq G V H L 1 (il). Note that, 
thanks to Remark l2.21 in this new regularity framework, (|3. 41) is nothing else than J*(i?o,n) < J*(^o), 
which as been proved above. 

Thus, we only need to prove the /^-convergence. To this end, we have to be a bit more careful. 
First, we define the Banach space I-fnL 1 )!!), endowed with the norm 

II • Ik : = II • Ik' + II • |U*(n), 

and introduce the unbounded linear operator A on X defined as Av :— Av with domain 

D{A) := {v G V : Aw G L 1 (0)}, 

where A is the usual distributional Laplace operator. Then, we have that A is an accretive operator 
on the space X. Indeed, by [H Prop. II. 3.1], this corresponds to checking that, if A > and 

Xi+XAxi=fi, i = l,2, (3.10) 

for fi G X , then 

\\xi-x 2 \\ x < ||/i-/a||jf. (3.11) 
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Actually, the analogue of p. lip w.r.t. the V^'-norm can be obtained by testing the difference 



xi - x 2 + XA( Xl - x 2 ) =fi-h (3-12) 

by (/ + A)~ 1 (x\ — X2) where I is the identity mapping of H (and, hence, / + A : V — > V' is the 
Riesz isomorphism). On the other hand, the ^-analogue of (|3.1ip is obtained by testing (|3.12l) by 
signal — X2) and applying the Brezis-Strauss theorem [9j Lemma 2]. Moreover, we have that D(A) 
is dense in V n L 1 ^). To see this, let us take z <eV (~) L 1 ^). Then, setting 

p n (x - y)z(y) dy, (3.13) 



n 

where p n is the standard mollifer, it is clear that z n is smooth (hence, in particular, it belongs to 
D{A)). Moreover, the convergence z n — > z in L 1 (f2) follows from standard properties of convolutions, 
while the convergence z n — >• z in V' follows from the density of H in V and from the fact that the 
mapping z 1— > z n is a contraction w.r.t. the V'-noim (this may be verified for z £ H by using Fubini's 
theorem and then extended to V by density). These facts permit to apply [U Prop. 3.2 (e)], which 
gives exactly the convergence property (|3.3[) . which concludes the proof. I 

Thus, taking $o,n as an initial datum for equation (|2.3[) (while the initial datum Xo is kept fixed), 
existence of a corresponding solution (d,u : X) is guaranteed by Theorem 12.11 Our aim will be now 
that of removing the approximation of the initial datum letting n /*• +00. With this aim, we start 
by recalling a couple of basic a-priori estimates. The procedure is detailed only in the conserved case, 
the differences occurring in the non-conserved case being pointed out at the end. For the meanwhile, 
we will not emphasize the dependence on n in the notation. 

Energy estimate. We test (|2.3[) by 1 + u, (|2.6[) by w and (|2.7p by X t . This formal procedure will 
be justified at the end, when we prove (|2.3ip . We obtain 

-^+||V M || 2 + ||VH| 2 =0, (3.14) 
at 

where £ was defined in (|2.20l) . Using also the properties of A, we immediately get 

\\£ Hl°°(o,t) + l|Vu|| L 2 (0T . ff) + ||Vw|| L 2 (0T;ff) + ||Xi|| L 2(o,T;V') < cE . (3.15) 

Here and below, the letters c and k will denote generic positive constants, independent of initial data 
and of time, whose value possibly varies on occurrence, k being used in estimates from below. In 
particular, the above estimate is uniform with respect to T. 

A generalized Poincare inequality. To estimate the full F-norm of u (and not just the iJ-norm 
of its gradient), we need a proper form of Poincare's inequality (cf., e.g., |19J Lemma 5.1] for a similar 
tool), which we prove just for the sake of completeness: 

Lemma 3.2. Assume £1 is a bounded open subset ofR d . Suppose v 6 W 1 ' 1 ^) and v > a.e. in ft. 
Then, setting K :— J n (\ogv) + , the following estimate holds: 

\\ v \\ Ll{n) <\n\e c ^ K + ^\\Vv\\ Llm , (3.16) 

the constants C\ and C2 depending only on Q. 

Proof. First of all, we recall that for any function z € W 1 ' 1 ^) such that |.Eq| > (with E a := 
z(x) = 0}) the following Poincare type inequality (see [3TJ Lemma 5.1, pag. 89]) holds: 

||«IUi<n) < Tp-rllVslUip,), (3.17) 

where the constant C can be explicitely computed and depends only on 17. Now, let v be a function in 
the hypothesis of the Lemma. Set K ;= / n (log v) + and note that, thanks to the Chebychev inequality, 
we have, for any fixed N > 0, 

\{xeft:(logv) + >N}\<^, (3.18) 
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and consequently 

\{ x€ Q: V> N}\< ] ^. (3.19) 
Thanks to ([3~18| and (l3~T9)l . we can fix N = N(K, ft) = e^T in such a way that 

| {a; e : v < N} \ = |fi| - | {x e Ct : v > N} | > (3.20) 
As a consequence, the inequality (|3.17p . with z = (v — N) + and (I3.20[) . entails 

|| W || Ll(n) < iV|0| + \\(v - iV) + || L1( o) < m\ + <jv}| l|V(u " ^) + H^(n) 

< |n|eM + 2C|n|||V«||£ 1{ n), (3.21) 

which is (|3TTC)1 with C x = 2/|Q| and C 2 = 2C, C being the constant in (|3~T7)) . I 

Consequences of the energy estimate. Using the above lemma, (|3.15[) additionally gives 

ll«||x»(t,t+i ; v) < Q(Eo) Vi >°- ( 3 - 22 ) 
Applying standard techniques to system (|2.6p - (|2.7[) . we also have 

\\ x \\L 2 (t,t+i;H 2 {n)) < QQEo). (3.23) 

Moreover, testing (|2.7p by X — Xq and proceeding, e.g., as in the Appendix of [23], it is not difficult 
to arrive at 

\\b(X)\\ L 2^ t+1 , H} < Q(E ), (3.24) 
whence a comparison of terms in (|2.7p and estimate p. 151) also give 

IMMm+W < Q(Eo). (3-25) 



Passage to the limit. We will now let n /• +oo, still referring to the conserved case. With this 
aim, we rename as (i9„,u„,X n ) the solution to the n- approximation. By p. 151) and (|3.22[) - p.25[) . we 
then have, for any T > 0, 

u n -> u weakly in L 2 (0, T; V), (3.26) 

X n ^X weakly in iJ 1 (0,T;U')n J L oo (0,T;U)n J L 2 (0,T;iJ 2 (fi)), (3.27) 

w n ^w weakly in L 2 (0,T;V), (3.28) 

b(X n )^b weakly in H X {Q,T;H), (3.29) 

for suitable limit functions u, X, b. Then, the Aubin-Lions compactness Lemma and the usual mono- 
tonicity argument [4j Prop. 1.1, p. 42] permit to see that b — b(X) a.e. in (0, T) x O. Moreover, the 
above relation suffice to pass to the limit in system (12.6[) - (|2.7p . 

Taking the limit in (|2.3|) and in (|2.4p is a bit more involved. Actually, (|3.15p and a comparison 
of terms in (|2.3p give 

K,t -> ^ weakly in i 2 (0, T; V'). (3.30) 
Then, integrating in time and using the ^'-convergence in (|3.2|l . we obtain more precisely 

i?„->tf weakly in if 1 (0,T;V). (3.31) 

This is sufficient to take the limit of equation (|2.3p . but not of (|2.4p . Actually, to identify $ in terms 
of u, we have to work a little bit more. Namely, we have to integrate (12.31) with respect to time both 
at the n-level and in the limit and then test, respectively, by u n and by u. Notice that, even at the 
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limit level, the use of u as a test function is guaranteed by the fact that u G L 2 (0, T; V) and all terms 
in flED lie at least in L 2 (0,T; V). 

Then, at the n-level, we obtain 

(i?„,w„) + -||V(l*u n )(i)|| 2 = / (tio,„ + X ,„-X n ,u n }, (3.32) 
z Jo 

where * denotes convolution in time. Taking the supremum limit in the above relation and comparing 
the result with the limit equation, we obtain 

limsup / (■&,!, u n ) < / (3.33) 
n/'+oo Jo Jo 

Since 

i? n = G 9^pJ(tiJ a.e. in (0, +oo), (3.34) 

Un 

(notice that we used here property (|2.16p ). relations (|3.26p and (|3.3ip and the standard monotonicity 
argument [H Prop. 1.1, p. 42], applied here in the duality pairing between V and V, permit to obtain 
(|2.30l) . which concludes the proof of existence in the conserved case. 

Differences occurring in the non-conserved case. At the level of estimates, the only relevant 
difference is in the energy relation, which is now obtained testing (|2.3[) by 1 + u and (|2.5D by Xt . This 
gives 

^+||V W || 2 + ||X 4 || 2 = 0. (3.35) 

Thus, we have the f/-norm of Xt rather than the T^'-norm on the left hand side (and consequently we 
obtain ([2^27)1 in place of ([2~56|) ). Estimates ((3~22|) - (|3T24| hold without variations while (j33i5l) makes 
no longer sense. Notice that p.24[) can now be obtained testing directly (|2.5[) by b(X). The passage 
to the limit is analogous. 

Proof of (|2.31[) . We first observe that (|2.30[) is equivalent to 

uedyyJ*^), (3.36) 

almost everywhere in (0,T). Then, the standard integration by parts formula [3 p. 73], applied in 
the duality between V and V, gives 

J*( 1 ?)6 4C([0 ) T]), <0 t ,u) = ^.r(0). (3.37) 

Thanks to this formula, in the non-conserved case for any ^'-solution we are allowed to test (12.31) by u 
and (I2.5P by X t . Integrating over (0,t) for arbitrary t > 0, we obtain exactly (|2.31[) . In the conserved 
case, instead, we have to test (|2.3p by u, (|2.6I) by w and (|2.7I) by X t . Note that this is still possible 
for any ^'-solutions. Indeed, thanks to (12.27)) and the properties of A, we have that w G L 2 (0, T; V). 
Thus, w can be used as a test function in (I2.6P (which is a relation in L 2 (0,T;V')) and Xt can be 
used as a test function in (|2.7I) (which is a relation in L 2 (0,T;V) thanks to the above discussion). 
However, we have to notice that, while 

AX + b(X) G L 2 (0, T; V), (3.38) 

it is not expected to be true that, separately, AX G L 2 (0,T; V) and 6(X) G L 2 (0,T; V). Nevertheless, 
as shown, e.g., in [29l Lemma 4.1], property (13.38[) is sufficient to prove that 

(x u AX + b{X)) = (I|VX| 2 + 6(X)) (3.39) 

almost everywhere in (0,T). Thus, we still have (12 .3 If) . of course with ||Vw|| 2 in place as ||X t || 2 . 

Proof of uniqueness. It works exactly as in the standard case (so, we just sketch it for the conserved 
model). Namely, we can take a couple of solutions Ui, Xi), ($2) 1*2, X2) starting from the same 
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initial datum, write the system for both solutions and take the difference. Then, setting (i?, u, X) := 
($1, Ui, Xi) — ($2, U2, X2), we integrate (the difference of) (|2.3[) in time and test it by u. Moreover, we 
test (the difference of) (|2.6p by A~ 1 X (note that X has zero-mean value, so A -1 is well-defined) and 
the difference of (|2.7[) by X. Collecting everything and noting that two couples of terms cancel, we 
obtain 

<0,u) + ~(||V(1 * u)\\ 2 + \\X\\ V ,) + ||VXf < ||X|| 2 < i||VX|| 2 + c\\Xf v , (3.40) 



Noting that, by monotonicity, (§, u) — — 1)2, u\ — U2) > 0, the thesis follows then from Gronwall's 
lemma. 

3.2 Proof of Theorem 12^41 

We start by deducing an additional a-priori estimate. As before, we present it just in the conserved 
case, the variations in the non-conserved case being given at the end. 

Second estimate — local version. We test (|2.3[) by tut = tdt/'d 2 and add the result to (I2.6[) 
multiplied by twt- Then, we add also the time derivative of (|2.7[) multiplied by tXt- We obtain 



^(^l|v.f + ^||vHI 2 ) + t/J + t||vx t f + ^ 2 ^(x)x 



Then, noting that 



<i(||V U || 2 + ||VH| 2 )+^H 2 - (3-41) 

t\\X t \\ 2 < ^||VX t || 2 + c<||VH| 2 , (3.42) 
integrating f|3.41[) between and t 6 (0, 1] and taking advantage of (|3.15[) . we obtain 

t\\Vu(t)\\ 2 + t\\Vw(t)f +t\\X t (t)\\ 2 v , + f S ||VX t ( S )|| 2 ds<cE . (3.43) 

Jo 

Second estimate — global version. We test (|2.3p by u t = i9t/$ 2 and add the result to (I2.6[) 

multiplied by w t . Then, we add also the time derivative of (12.71) multiplied by X t . Proceeding as 
above, we obtain 

32 



^(^l|Vu|| 2 + i||V W || 2 )+ f ^ + i||VX t || 2 + / fe'(X)X 2 <c||V W || 2 . (3.44) 
Integrating between 1 and t > 1 and recalling (13. and (|3.43p , we infer 

/oo 
||VX t (i)|| 2 dt<Q(E ). (3.45) 

Using again the logarithmic Poincare inequality (Lemma 13. 16[) , we also have 

IMIl~(i,oo;V0 < Q(Eo), (3.46) 

i.e., the first (1232)) . 

In the non-conserved case, the procedure is similar. In place of (|3.45[) we rather obtain 



||Vu|| i0 o (1>oo;J?) + ||X t || 1 «. (li00 . JI) + J ||VX t (i)|| 2 dt<Q(E ). (3.47) 

As a further consequence, we can look at equation (|2.7I) in the conserved case f (|2.5p in the non- 
conserved case, respectively). Thanks to estimates (I3.45[) - (|3.46[) for u and w (respectively, to estimate 
(|3.47p for u), applying standard regularity results for elliptic equations with monotone nonlinearities, 
we then obtain (|2.33l) . 

Asymptotic uniform regularity of u. Our aim is now to show the second (|2.32p . The key step is 
represented by the following lemma: 
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Lemma 3.3. Let u be a solution of the problem 

ti t +Au = f, § = -l/u (3.48) 
over the time interval (S, S + 2), where we additionally assume that 

H M llL3(S,S+2;L 3 /2(n)) < M, ||/||l2( S ,s+2;L 3 + = (0)) < F, (3.49) 

for some (given) constants M > 0, F > and some e > 0. Moreover, let us assume that 

u(S) e L 1 ^). (3.50) 

Then, we have 

IMU-((s+i,s+2)xn) < Q(F,M, ||«(5)||i). (3.51) 

Proof. We test l|3.48|) by — |w| p+1 , where p > 1 will be specified later (although needs 
not necessarily be an admissible test function, the procedure could be easily justified by truncation 
arguments, we omit the details). This gives 

1 d MI5+ ^y^HVN^H 3 < / |/||.r. (3.52) 



pdt" " p (p + 2) 

We then set r := to be the conjugate exponent of 3 + e. Then, multiplying by p, we can estimate 
the right hand side as 

P+m =«llflL, AUMPf 1 



P / |/lhr 1 <p||/]] 3+£ |||< + 1 r =p||/||3 +£ ||«||^ 1 +1) . (3.53) 
Then, in order to recover the full y-norm from the gradient term, we add 

IIm^IIHnS! (3-54) 

2 

to both hands sides of (I3.52p . Integrating (I3.52p over (r, t), for t a generic point in (t,S + 2) and 
choosing, for the first iteration, p = 1 and r = S, we obtain 

W U W P L^(T,S+2:LP(fl)) + IMIi>+ 2 (T,S+2;L 3 P+ 6 (n)) 

rS+2 

< 4 u (t)\\ P lp(q) +cp||/||L2 (SiS+2;L 3+ e(a)) ||M||^t( 1 P+ i) (r . 5+2;i .(p + i) (n)) +c / ||u(s)|fgf ds 

Jt 2 

< c P F \\ u \\ P L2l + i ){ ^ s+2 . L r {P+ i ){n)) + Q{M, \\u(S)\\ p ), (3.55) 

where in the last inequality we took advantage of p.49[) using that p = 1 and t = S. 

Being non-restrictive to assume that u > 1 almost everywhere (otherwise, we can replace u 
with max{u, 1}), we can then define 

Jp '■— IMIl~(t p ,S+2;Lp(Q)) + IMIlp+ 2 (t p , S+2,L 3 J>+8(n)) ' (3.56) 

where, for now, we take t v = r = 0. Then, by interpolation we obtain 

l|w||i2(»+i)(T 1 ,,S+2;i'-<«+l)(n)) < H U llL°°(r p ,S+2:LP(n))ll U llLP+2(T p ,S+2,L 3 P+6(n))' ( 3 - 57 ) 

for some a G (0,1)- Then, raising to the power p and using the Young inequality with exponents 
P = 1/ot e Q = 1/(1 — a), we get 

ll U llL 2 (i+ 1 )(r p ,S+2;L''(<!+i)(n)) - a IHlL° c (T p ,S+2;LJ>(fi)) + - °0 H U lliP+ 2 (T p ,5+2,L 3 P+ 6 (n)) ' (3.58) 

which implies, upon dividing by max{a, I — a} (that is different from and 1) 

Jp > IMlL2(,+i)( Tpj 5 + 2 ; L>'(<!+i)(fi))> (3.59) 
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where the index q and the interpolation exponent p are given by the system 



p+2 2(5+1) ' 

£ _i_ 1 ~P _ 1 
. p 3(p+2) r(q+l) ' 



Dividing the second equation in (|3.60[) by the first one, we actually have 

'p l-p \P + 2 2 



p 3(p + 2)/l-p r 



(3.60) 



(3.61) 



whence 

2E±2_H ' *„ (3.62) 
p 1 — p r 3 

and it is easy to compute 

K E = Jiif . (3.63) 
9 + 3e v y 

From (|3.62p and the first equation in (|3. 601) . we then have 

1 + ^+2^ 



Being 



we then obtain fron the hrst (|3.60p 



1 ~ P= -i T p F ' ( 3 ' 65 ) 



1 p + 2 1 / p 

- 1 =o( 1 + ^^ (P + 2)-l 



21-/9 2V p + 2 

^s + 1 9 + 5e , . 

= ~^~ P = oTs? =: HP > (3 ' 66) 

where, obviously, if = H (e) > 1 whenever e > 0. 

Given that po = 1, let us set, by induction, Pi+i = Hpi = H %+ , Then, let z > and let us 
rewrite p.55[) by taking p = Pi+i and r = Tj+i (the latter will be chosen below). Setting also, for 
brevity, Ji := J Pi , we then obtain, thanks also to (|3.59p . 



J%? < c\\u(n +1 )\\^ i+ cp l+1 FJ?^ +1 +c f \\u{s)\\^ + \%l ds. 
Now, let (for instance), for i > 1, 

6 °° 
Oi = o~^t ; so that > Ui = 1. (3.68) 

4 — 1 

Then, we observe that, given Tj, we can choose G (Tj, + o"j+i) such that 

||«(rn-i)ll?J +I < ~r~ r +ffl+1 IN*)II£* +1 df 

< cz 2 £ +aZ+1 \\u(t)\\ p 3 l pi+2) At < d?Jf , (3.69) 

where we used that Pi+\ < 3(pi + 2). 
Analogously, we have that 

5+2 \HsmZ\tl ds < 2c\\u\\ p ^ s+2 . LPi{n)) < cJp +1+2 . (3.70) 

Tj+1 2 



(3.67) 
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Collecting (P^ - ljXTD)) . (j3~57) gives 



J&F < ci 2H JT +1 + cp l+1 FJ^ 1+1 + cJf +1+2 . (3.71) 

Thus, we finally obtain 

JUT < c{i 2H + H i+1 F + 1) jf +1+2 . (3.72) 
Thus, setting ry 4 := (H* + 2) /if*, we have 

J+i < Sf _<I+1) Jf +1 , where fl, := c(z 2ff + iP +1 F + l) (3.73) 

whence a simple induction argument (see, e.g., [30]) permits to obtain (|3.5ip . I 

Conclusion of the proof. To obtain the uniform boundedness of u it is now sufficient to notice that, 
by the first (|2T3"2"|) . u\t) G L 6 (fi) for all t > 1. Then, (|3~STJ| holds. Moreover, (|3T35|) are a consequence 
of (|3.45p (or, in the non-conserved case, (|3.47[1 ). which gives the required regularity for / = — X t , and 
of (|3.46|) . which gives the required regularity for u. The second (|2.32[1 is then a consequence of the 
lemma (applied with the choice of 5 = 2). It is also worth noting that, at the level of ^'-solution, the 
identification property in (|3.48[) needs not hold in the strong (pointwise) form (but just in the "weak" 
sense (|2.30p ). However, one can apply Lemma [3.31 at the ro-regularized level, and then pass to the 
limit noting that the procedure yields n-uniform estimates. 

To complete the proof we have to show (|2.34|) in the non-conserved case. Of course, such a 
property is significant only in the case when /, the domain of 6, does not coincide with the real line, 
i.e., we are in presence of a singular potential (like the logarithmic one (jl.ip ). Otherwise, (|2.34|) is 
(also in the conserved case, of course) an immediate consequence of (|2.33p . 

That said, let us prove the upper bound, the lower one working in a similar way. Being, by 
the second (|2.32p . \u(t, x)\ < U for some U = Q(Eo) and a.e. (t, x) € (2, oo) x f2, we can then apply 
the comparison principle to (|2.5p . This gives that X is bounded from above by the solution X + to the 
forward Cauchy problem 

X+ + b(X+) - X + = U, X+(2)=i+, where 1+ := sup/. (3.74) 

Actually, by the last (|2.9p . lim T - / / sup / b(r) ~r — U = +oo. Thus, X + (t) < I + — 5 for all i > 3 and some 
S > 0. Moreover, this bound is uniform in time. Then, (|2.34p is a consequence of the comparison 
principle. This concludes the proof of the theorem. 



3.3 Proof of Theorem 12751 

We start by giving the proof in the non-conserved case, the variations occurring in the conserved case 
being outlined at the end. Then, we know that any ^'-solution (i?, u,X) satisfies the energy equality 
(|2.3ip . Analogously, if (t? n ,u n ,X n ) is the approximating solution constructed in the existence proof, 
the analogue of (|2.3ip reads 

/ (||V Un || 2 + ||X„, t || 2 ) = J*(tfo,«) - J*(4 n (t)) + (0o,», 1) - <!?»(*), 1) 
Jo 

+ 1 (^o + ^|VX | 2 +S(Xo)-ix 2 ) - f (x n {t) + \\VX n {t)\ 2 +b(X n (t)) - i*„(i) 2 ). (3.75) 

Our task is now to compute the supremum limit of (|3.75l) and compare it with (|2.3ip . Then, we firstly 
observe that, by (|3.27p . the Aubin-Lions lemma, and lower semicontinuity of &, 

/ + hvx(t)\ 2 +bWt)) - ^X(t) 2 ) < Hminf / (x n (t) + \\VX n (t)\ 2 +b(X n (t)) - \x n {tf). 

J Q V 2 2 J nSoa J n \ 2 2 J 

(3.76) 

Analogously, using convexity and lower semicontinuity of the functional J* w.r.t. the y'-norm, (|3.2I) . 
and the fact that, by (|3.31l) . "& n [t) tends to weakly in V for all t € [0, T], we obtain 

J*(tf(i)) - (0 , 1) + W*), 1) < Hminf ( J*(i?„(t)) - (^ ,„, 1} + (l? B (t), 1)). (3.77) 
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So, it remains to prove that 

limsup J*(i? ,n) < J*("&o)- (3.78) 

n/oc 

The proof of this fact is actually a bit more involved. We prepare a Lemma 

Lemma 3.4. Let v G V H L 1 ^) such that j*(v) G L^O), j and j* being given by (|2J2|l . Then, 



■/» - / j»- (3-79) 

Proof. Let 2 belong to the domain of J, namely let z G with J(z) < +00. Then, for a.c. x £ f2, 
by definition of subdifferential in R, we have 

^(i)«(i)<j(«(i))+j'(»(i)). (3-80) 
Integrating over fi, we would formally get 

zjj < J(z) + / j*(v). (3.81) 



However, the integral on the left hand side could make no sense since the function zv could not belong 
to T 1 (i7). Nevertheless, it is simple (see, e.g., [51 Lemma 2.2] and [H Lemma 2.1]) to see that, in place 
of (ETgTj) . there holds 

(z,v) <J(z)+ f j*(v). (3.82) 
Jn 

Passing to the supremum w.r.t. z varying in the domain of J, we then get the < sign in p.79[) . 
To prove the converse, we first let, for v as in the statement and N G (1, 00), 

z N := -max{iV _1 ,min{w _1 ,iV}}. (3.83) 

Then, for e G (0, 1), we regularize by singular perturbation as in (|3.1[) : 

z e ,N + cAz CtN = z N . (3.84) 

Then. 



J*(v) > (v,Z €j n) - J(Ze,N) = / VZc.N - J{z<l,n) = I VZ e<N + / log(-Z e ,Jv) 

Jn Jn Jn 

> / vz,. N + / log(-^jv) ^ eN ° / vz N + / log(-Zjv) 
Jn Jn Jn Jn 

f (-l + \ogv- 1 )= f j», (3.85) 
Jn Jn 

where the first inequality follows from definition of conjugate function, the second equality from the 
fact that z €i n is smooth and bounded, the third equality is trivial, the fourth inequality comes from 
(|3.84[) . the convergence e \ from standard properties of elliptic systems, and the convergence 
N /~ 00 from Lebesgue's theorem. This proves the > in (|3.79p and the lemma. I 

Then, by definition of conjugate function (recall (|2.12p ). 

J*(fio,n) = SUp (#0,7i, «) ~J(v)< SUp (l?0,n,w) - J(v) 
v£V vEH 

= [ -1 - logi? ,n < [ -1 - logi?o, (3.86) 
Jn Jn 

where the last inequality follows from (13.41) . In particular, we have that —1 — log$o G i 1 (fi). Thus, 
applying the above Lemma, we obtain 

J*(flo) = / (-1-logtfo), (3-87) 
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whence, computing the supremum limit of p.86[) . (|3.78|l follows. Thus, we finally end up with 



limsup f (||Vit„[[ 2 + ||X r 

n/oc Jo 



< 



iVull 



ll*t| 



(3.88) 



whence, recalling (I3.26[) - (|3.29[) . we get in particular 

Vu„, Xn,t Vu, X t strongly in L 2 (0,T; H). 

As a final step of our procedure, we shall prove that is a Cauchy sequence in C°([0, T]; 

Actually, writing equation (|2.3p for a couple of indexes n and m and taking the difference, we obtain 



(3.89) 



$n,t ~ $m,t + A(u n — u rn ) — X n ,t — X m> t. 



(3.90) 



Thus, testing by sign(t9„ — z9 m ), noticing that, by monotonicity, sign($„ — d m ) = sign(w„ — u m ), and 
applying the Brezis-Strauss theorem [9j Lemma 2], we arrive at 



^ll^n ~ ^m||Li(0) < \\X n ,t - X m ,t||£i(n), 



(3.91) 

whence, integrating in time and using the strong convergences p. 21) and (|3.89p . we end up with 

0„-H? strongly in C°([0, T]; i x (0)). (3.92) 

Thus, in particular, we have that, for all t E [0, T], E i 1 (0) D <9y,y J(u(t)), whence the pointwise 
identification (12 .4[) follows from (|2.17l) . This concludes the proof in the non-conserved case. 
Conserved case. To conclude the proof, we outline the differences occurring in the conserved case, 
which only regard the above i^-argument. Actually, the convergence of X n t in f|3 . 89[) is now replaced 
by 

Vw„ -> Vw, strongly in L 2 (0,T;H). (3.93) 

Of course, thanks to the properties of A this also gives 

X n<t -> X t strongly in L 2 (0,T; V'), (3.94) 

which, however, is not sufficient to proceed as before. On the other hand, we can rely on estimates 
(jiRBl and (|3~4"Tj) which tell us that 



t^X 



n > t \\L 2 (0,T;V) 



4- Wt 1 / 2 ^ t II < r 

-r \\ L a ^A\l^(o,t-v) - c - 



Thus, by interpolation, 



k^x 



n til -J— < C 



(3.95) 

(3.96) 



for all e E (0, 1). Coming back to (|3.91l) . we now have that, for all t E [0, T], 

l-e 1 

\\X n ,t — ^m,t||i 1 (fi) ^ c||X„ it — Xm : t|| < c||X„ :t — X m> t||y, E \\X n j — X my t\\ H i 



'(H) 



< ct 2 < 2 -> ||X„, t - X m , t ||fT • i 2 < 2 -> ||X B)t - X ro , t ||^T-. (n) 
Then, integrating over (0,T), we arrive at 

||X„, t -X 



l -m,t||L 1 ((0,T)xa) 



< c 



< c 



* 2(2 - E) \\Xn,t ~X m ,t 
T 



L3(0J) 



t^T X„, t -X m , t 



(3.97) 



L 4 (0,T) 



2 «(!-«) 

r^)||x„, t -x m , t ||^ 2 - £) 



ii 



< c i ^f 5 ^ 



L 2 (0,T) 



X 



m,t\\yi 



4(1-0 
3(2-0 



L 2 (0,T) 



c||X nj i X mj i|| 8(l-e) 5 

i 3(2-.) (0,T;V) 



(3.98) 



where the second inequality is a consequence of (|3.96p . the thirds follows from Holder's inequality, 
and the fourth holds provided that we take e so small that 4 < 3(2 — e). In particular, using (|3.94p 
(note that 8(1 — e)/3(2 — e) is smaller than 2 for e as above), we obtain that the right hand side tends 
to for large m and n. At this point the proof goes on like in the non-conserved case. 
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3.4 Proof of Theorem 12771 

To start, we need to prove some further a-priori estimates holding under the additional assump- 
tion (j2.37[) . In particular, the key step will be that of showing that the L 3+£ regularity of $ is 
conserved uniformly in time. To show this, we will use L p -techniques in equation (|2.3[) (actually, with 
p = 3 + e). However, doing this will require some care since, due to the low regularity of initial data, 
the "forcing term" X t needs not belong to L 3+e (0) for small values of the time variable. However, 
we will see that the L 3+e -norm of X t (t) explodes, as t \ 0, in a way which is sufficiently slow for our 
purpose. As before, the proof is detailed just in the conserved case. That said, we start with the 

Third estimate — local version. To start, we test (|2.3I) by d 2+£ , to obtain 



-I 
elf 



whence, clearly, 

Thus, using that 
and that 

we obtain 



\l + + l<c\ \X^ 2+e <c\\X4, +s \\nl\l, (8.!.!)) 

^||tf||3 +e <c||X t || 3+e . (3.100) 

#§$i(fi) c£ 3+e (0) (3.101) 

— =aXl + (l-o)x(-l). with a = JL±ii (3.102) 

6 + 2e v ' y ' 12 + 4e K J 



A 9+5e 3-e 

- IL9ll 3 +s <c||X t ||^||X t ||^ 



dt 



<Q(E ){t^\\X t \\^)t-* 

<Q(E )(t\\X t \\ 2 v + t-^), (3.103) 

where (|3.43[) has also been exploited. Note that, for e G (0, 1) the latter exponent lies in (—1, 0). Notice 
also that in the non-conserved case the exponents are even better since it is sufficient to interpolate 
between V and H (rather than between V and V). Thus, integrating (13. 103|) between and 1, and 
using once more (|3. 431) . we infer 

l|0|U~(o,i;z3+. ( n)) < Q(Eo, ||^o||3+ e ). (3.104) 

Third estimate — global version. As before, we test (12.3[) by i} 2+£ . Taking now care also of the 
gradient term, we get 

d_ 

dt 



: ._ : -r«||W = || <c / \X t \d 2+£ . (3.105) 



Adding also the inequality (which is true thanks to (|3.15|0 

4^ \\l < c(l + < Q(Eo), (3.106) 



we then get 



^]]^ + K ||^||^<c||X t ||6||^|| 3/2 ||^|| 6 + Q(Eo) 



<c\\X t \\ 2 v \\^\\l /2 + ^p^\\l + Q(E, 
< c\\X t \\ 2 v p\\l + + l + ^||^||6 + Q(Eq) 
< C ||X t ||^||,?||^ + ^||^||^ + Q(E ). (3.107) 
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Thus, setting Y(t) := ||i9(t)||§ +e , we obtain from p.!07|) 

Y' + K<c\\X t \\ 2 v Y + QY- 1 i £ , (3.108) 
where we wrote Q in place of Q(Eo). Now, let us set 

:=max{r(t), (^^J, (3.109) 

so that it is clear that Z satisfies, 

Z'(t)<c\\X t (t)\\ 2 v Z(t), fora.e. t>l, and Z(l) < Q(E , \\MU+e) ■ (3.110) 
Hence, integrating between 1 and a generic t > 1 and recalling p.45[) . we arrive at 

||0|U-(l,oo;£»+-(n)) < Q(Eo, H^olU+e), (3.111) 
i.e., the global analogue of (|3.104|> . 

Asymptotic uniform regularity of -d. The key step is represented by the following counterpart of 
Lemma 13.31 



Lemma 3.5. Let u be a solution of the problem 

ti t +Au = f, tf=-l/u, (3.112) 
over the time interval (5, 5 + 2), where we additionally assume that 

||tf|U~(S,S+2;£i(n)) < A/, ||/||^(S,S+2;i3+e(n)) < F, (3.113) 

for some (given) constants M > 0, F > and some e > 0. Moreover, let us assume that $(S) 6 
L 3+e (fi) for some e > 0. Then, 

||0||i~((5+i,s+2)xn) <Q(^M,P(5)|| 3+E ). (3.114) 
PROOF. We test (13.481) by where p > 3 will be specified later. This gives 

Setting as before r := |±| and multiplying by p, we have 

P / l/l^" 1 ^Pll/lla+ell^-lr =P||/|| 3+£ |KH^ 1 _ 1) (3.116) 
and, in order to recover the full V-nonn from the gradient term, we add 

\\^\\l = \^Cl { ^y ( 3 - 117 ) 

Integrating (|3.115p over (r, t), for t a generic point in (r, 5 + 2) and choosing, for the first iteration, 
p = 3 + e and r = 5, we then have the analogue of (|3 . 55[) : 

ll^llx,°°(T 1 S+2;£*(n)) + ll^lli>- 2 (T,S+2;L 3 P-6(n)) 

/■S+2 

< c|k(r)f iP(0) +cpF||^f L + 1 p _ 1)(TS+2;Lr(p _ 1)(0)) +c jf 

< cpF\m^_ lKrjt . Lr(p _ 1Ka)) + Q(M, Huolli), (3.118) 

Now, the iteration scheme goes through similarly as before. Actually, in place of (|3.60p . we get the 
system 

[p + 3(p-2) r(g-l) ' 
18 



whence one computes, exactly as before, 



9 + 7e 

9 + 3e y ' 



1P ~ 2 ' l = i(l+-M B )(P-2) + l 



and, finally, 



2 1 - p 2 V p- 2 

= = 9T3i P = (3 ' 121) 

which, exactly as before, is larger than one. Hence, the procedure continues as before, with small 
variations in the numerical values of the indices. Of course, the L 3+£ regularity of the initial datum 
is used since we need to take p = po = 3 + s at the first iteration (for smaller values of p we get no 
summability gain from the gradient term). I 

Conclusion of proof. Thanks to estimate (|3.111[) , $ satisfies the first (|3.113p for any S > 1 (where 
M is the quantity on the right hand side of (I3.11ip which is independent of S). Analogously, we have 
$(S) £ L 3+£ for (almost) all S > 1. Moreover, combining Q3.15P and (|3.45p . we have the second of 
(|3.113p . still with F independent of S. We then conclude applying the above Lemma over the generic 
interval (S, S + 2), with S > 1. 

Remark 3.6. Of course, with (|2.32l) and (12.38)) at our disposal, equation (|2.3p is both nonsingular 
and nondegenerate. Consequently, we can prove, with standard tools, further regularization properties 
of solutions. In the non-conserved case, thanks to (|2.34l) . also the (possibly) singular character of b is 
lost. Thus, the smoothness of the solution is limited only by the differentiability properties of b. For 
instance, if b £ C°°, then also the solution is infinitely diffcrentiable for strictly positive times. 



References 

[1] S. Agmon, A. Doughs, and L. Nirenberg, Estimates near the boundary for solutions of elliptic 
partial differential equations satisfying general boundary conditions. J., Comm. Pure Appl. Math., 
12 (1959), 623-727. 

[2] N.D. Alikakos, L p bounds of solutions of reaction-diffusion equations, Comm. Partial Differential 
Equations, 4 (1979), 827-868. 

[3] C. Baiocchi, Sul/e equazioni differenziali astratte lineari del primo e del secondo ordine negli spazi 
di Hilbert (Italian), Ann. Mat. Pura Appl. (4) 76 (1967), 233-304. 

[4] V. Barbu, "Nonlinear Semigroups and Differential Equations in Banach Spaces", Noordhoff, 
Leyden, 1976. 

[5] V. Barbu, P. Colli, G. Gilardi, and M. Grasselli, Existence, uniqueness, and longtime behavior 
for a nonlinear Volterra integrodifferential equation, Differential Integral Equations, 13 (2000), 
1233-1262. 

[6] M. Bonforte and J.L. Vazquez, Positivity, local smoothing, and Harnack inequalities for very fast 
diffusion equations, Adv. Math., 2 (2010), 529-578. 

[7] H. Brezis, "Operateurs Maximaux Monotones et Semi-groupes de Contractions dans les Espaces 
de Hilbert", North-Holland Math. Studies 5, North-Holland, Amsterdam, 1973. 

[8] H. Brezis, Integrates convexes dans les espaces de Sobolev, Israel J. Math., 13 (1972), 9-23. 

[9] H. Brezis and W.A. Strauss, Semi-linear second-order elliptic equations in L 1 , J. Math. Soc. 
Japan, 25 (1973), 565-590. 



19 



[10] P. Colli, G. Gilardi, E. Rocca, and G. Schimperna, On a Penrose-Fife phase-field model with non- 
homogeneous Neumann boundary conditions for the temperature, Differential Integral Equations, 
17 (2004), 511-534. 

[11] P. Colli and Ph. Laurengot, Weak solutions to the Penrose-Fife phase field model for a class of 
admissible heat flux laws, Phys. D, 111 (1998), 311-334. 

[12] P. Colli, Ph. Laurengot, and J. Sprekels, Global solution to the Penrose-Fife phase field model 
with special heat flux laws, Variations of domain and free-boundary problems in solid mechanics 
(Paris, 1997), 181-188, Solid Mech. Appl., 66, Kluwer Acad. Publ., Dordrecht, 1999. 

[13] A. Damlamian and N. Kenmochi, Evolution equations generated by subdifferentials in the dual 
space of if 1 (ft), Discrete Contin. Dynam. Systems, 5 (1999), 269-278. 

[14] E. Feireisl and G. Schimperna, Large time behavior of solutions to Penrose-Fife phase change 
models, Math. Methods Appl. Sci., 28 (2005), 2117-2132. 

[15] G. Gilardi and A. Marson, On a Penrose-Fife type system with Dirichlet boundary conditions for 
the temperature, Math. Methods Appl. Sci., 26 (2003), 1303-1325. 

[16] M. Grun-Rchomme, Caracterisation du sous-differentiel d'integrandes convexes dans les espaces 
de Sobolev (French), J. Math. Pures Appl. (9), 56 (1977), 149-156. 

[17] A. Ito and N. Kenmochi, Inertial set for a phase transition model of Penrose-Fife type, Adv. 
Math. Sci. Appl., 10 (2000), 353-374 (Correction: Adv. Math. Sci. Appl., 11 (2001), 481). 

[18] N. Kenmochi and M. Kubo, Weak solutions of nonlinear systems for non-isothermal phase tran- 
sitions, Adv. Math. Sci. Appl., 9 (1999), 499-521. 

[19] N. Kenmochi, Neumann problems for a class of nonlinear degenerate parabolic equations, Differ- 
ential Integral Equations, 3 (1990), 253-273. 

[20] N. Kenmochi and M. Niezgodka, Systems of nonlinear parabolic equations for phase change 
problems, Adv. Math. Sci. Appl., 3 (1993/94), 89-117. 

[21] O.A. Ladyzhenskaya, V.A. Solonnikov, and N.N. Ural'ceva, "Linear and Quasi-linear Equations 
of Parabolic Type", Transl. Math. Monogr., Amer. Math. Soc, Providence, RI, 1968. 

[22] Ph. Laurengot, Solutions to a Penrose-Fife model of phase-Geld type, J. Math. Anal. Appl., 185 
(1994), 262-274. 

[23] Ph. Laurengot, Weak solutions to a Penrose-Fife model with Fourier law for the temperature, 
J. Math. Anal. Appl., 219 (1998), 331-343. 

[24] A. Miranville and S. Zelik, Robust exponential attractors for Cahn-Hilliard type equations with 
singular potentials, Math. Methods Appl. Sci. 27 (2004), 545-582. 

[25] R.H. Nochetto and G. Savare, Nonlinear evolution governed by accretive operators in Banach 
spaces: error control and applications, Math. Models Methods Appl. Sci. 16 (2006), 439-477. 

[26] O. Penrose and P.C. Fife, Thermodynamically consistent models of phase-Geld type for the ki- 
netics of phase transitions, Phys. D, 43 (1990), 44-62. 

[27] O. Penrose and P.C. Fife, On the relation between the standard phase-Geld model and a "ther- 
modynamically consistent" phase-Geld model, Phys. D, 69 (1993), 107-113. 

[28] J. Priiss and M. Wilke, On conserved Penrose-Fife type models, arXiv:1002.0928, to appear in 
Progress in Nonlinear Differential Equations and Their Applications. 

[29] E. Rocca and G. Schimperna, Universal attractor for some singular phase transition systems, 
Phys. D, 192 (2004), 279-307. 



20 



[30] G. Schimperna, Global and exponential attractors for the Penrose-Fife system, Math. Models 
Methods Appl. Sci., 19 (2009), 969-991. 

[31] R.E. Showalter, "Monotone Operators in Banach space and Nonlinear Partial Differential Equa- 
tions" . Mathematical Surveys and Monographs, 49. American Mathematical Society, Providence, 
RI, 1997. 

[32] J. Simon, Compact sets in the space L p (0,T;B), Ann. Mat. Pura Appl. (4), 146 (1987), 65-96. 

[33] J. Sprckels and S. Zheng, Global smooth solutions to a thermodynamically consistent model of 
phase-held type in higher space dimensions, J. Math. Anal. Appl., 176 (1993), 200-223. 

[34] J.L. Vazquez, "Smoothing and decay estimates for nonlinear diffusion equations", Oxford Uni- 
versity Press, Oxford, 2006. 



21 



